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ABSTRACT 



A three-level quasi-geostrophlc atmospheric model on a 
mid-latitude B-plane with a single wave in the x-direction 
and a thick lower layer was developed. This model was used 
to determine the disturbance structure in y for linear inter- 
actions with a constant mean v^^ind and nonlinear interactions 
with a varying mean wind. Simple heating to maintain the 
mean wind and friction were included in the nonlinear case. 

The linear experiment produced exponential growth rates 
in the disturbance field that were independent of the lower 
layer thickness. These viaves were baroclinically unstable. 
The mean v;ind field showed large modifications for a Rossby 
deformation length of 75 KM at 104 Jupiter days. 
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I. INTRODUCTION 



The meridional belts (dark bands) and zones (light bands) 
of Jupiter have been recorded since this planet was first 
observed through a telescope. The features found in these 
regions are semi-permanent and have a variety of rotation 
rates. Peek (1958) assembled this information to shov/ the 
variation of zonal velocity with Jovian latitude. Evidence 
was found by Chapman (19^9) for anticyclonlc shear to be 
associated with the zones and for cyclonic shear to be 
associated with the belts, but the correlation was Imperfect. 

Ingersoll and Cuzzi (19^9) have investigated the baro- 
troplc Instability of the observed motions at mid-latitudes 
using data from Peek (1958). Geostrophic balance and 
symmetric temperature difference between the belts and zones 
were assumed. The validity of the assumptions was verified 
by the excellent agreement between predicted and observed 
zonal velocities. It was also shown that the relative 
vorticity gradient U" approaches, but does not exceed the 
planetary vorticity gradient 6. This is the necessary 
condition for barotrophic instability. Thus barotrophlc 
instability probably limits either the belt wave number 
or the amplitude of the zonal flow. It was further shown 
that the zones are warmer than the belts. 

Stone (1967) applied the baroclinic stability theory 
to the dynamics of the Jovian atmosphere. He tested the 
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hypothesis that the zonal motions in Jupiter’s atmosphere 
are thermal v;inds and that the latitudinal cloud bands are 
caused by baroclinic instabilities under non-geostrophic 
conditions. Stone (1966a) showed that at least three 
different kinds of baroclinic instability can occur when 
d6/dZ > 0. The kind of instability that dominates depends 
upon the value of the non-dimensional Richardson number 

SI - I 



where g is the acceleration of gravity, 0 is the potential 
1 9 0 

temperature, ^ ^ is a representative vertical stability, 

9 U 

and 7^7 a characteristic vertical shear of the wind. The 

o Zt 

Richardson number is a characteristic ratio of work done 
against gravitational stability to energy transferred from 
mean to turbulent motion. He assumed that the observed 
v;inds on Jupiter are due to symmetrical Instabilities with 
no longitudinal variations. This type of instability domin- 
ates in the range' 0.25 < Ri < 0.95* Allowance for a deep 
atmosphere was made by not using the Boussinesq approximation 
for compressibility. It was shown that in a deep atmosphere 
symmetric instabilities form in the specified range of Ri 
and that these instabilities may extend over many scale 
heights. Thus the zonal winds and the latitudinal cloud 
bands may be explained qualitatively by assuming the presence 
of north-sourth temperature gradients caused by differential 
solar heating. These gradients will drive thermal winds in 
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the zonal direction, and under appropriate conditions, 
symmetric instabilities can grow, generating latitudinal 
cloud bands. 

Gierasch, Ingersoll, and Williams (1973) demonstrated 
that a cloudy planetary atmosphere at rest is unstable to 
disturbances of large horizontal scale. The energy source 
for this instability is the change in radiative heat flux 
associated v;ith a vertical displacement near the emitting 
level. The important feature of the model v;as that the net 
heating rate Q is proportional to the vertical displacement 
6Z of the cloud from its equilibrium position. Quasi-steady 
geostrophic motions or short period oscillations, representing 
inertia-gravity waves may become unstable in the presence 
of radiative heating. Gierasch and Goady (1969) found that 
the radiative time constants of Jupiter are longer than the 
dynamical time constants in the earth's atmosphere; there- 
fore, the geostrophic mode was applied, since the instability 
can be quasi-geostrophic . Gierasch, Ingersoll, and Williams 
( 1973 ) derived a value of 19,000 KM for the y-scale wavelength 
at 23 ° latitude which corresponds to the limiting amplitude 
for an axlsymmetric geostrophic instability where the wave- 
number in the x-direction k =0. This y-length is essen- 
tially equal to the observed belt-zone pair at this latitude. 
They also computed the value of the radiative time constant 
For an isothermal atmosphere X^ ^ was approximately 
equal to l50 years. It v;as also observed that even if the 
lapse rate differed from the adiabatic lapse rate by only 
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1 % , the radiative time constant would be about 1.6 years. 
These results are consistent v;ith the observed lifetimes 
of atmospheric features on Jupiter. 

Gierasch (1973) studied the radiative instability of 
a cloud deck on an equatorial 3-plane as discussed by 
Gierasch, Ingersoll, and Williams (1973). The model was 
based on the assumption that a precipitating condensable 
atmospheric constituent can be concentrated in the atmos- 
phere by vertical motion. If the condensing constituent is 
radlatively Important, it can cause radiative warming at 
the locations v;here it is concentrated. As a result, the 
heating is related to vertical motion and reinforcement is 
possible. They considered two values of the Rossby 
deformation length 



where g' is reduced gravity, f is the Coriolis parameter, 
and H is a characteristic height. Static stability esti- 
mates by Levels (1969) and Stone (1972) gave characteristic 
lengths of 785 KM and 78.5 KM, respectively. It was con- 
cluded that the radiative Instability mechanism driven by 
the dynamical concentration of a condensable constituent, 
produces instabilities with the wavelength and zonal symmetry 
of Jupiter's bands. Nonsymmetric instabilities with a 
shorter wavelength also exist. Finally, the structure of 




( 1 . 2 ) 
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the modes is in agreement with Ingersoll and Cuzzi's (1969) 
conclusion that cloudy zones are v;armer than belts. 

A three-level quasi-geostrophic model Incorporating a 
deep lower layer, as suggested by Gierasch (1973), was 
developed based on the Charney and Phillips (1953) multi- 
level model. The model was modified by Implementing Lorenz's 
(i 960 ) technique of replacing the dependent variable with 
Fourier coefficients. 

The objective of this study was to define the linear 
disturbance structure in y and to see if the belt-zone 
wind growth rate due to quasi-geostrophic barocllnic 
instability decreases with greater lov;er layer depth. Non- 
linear interactions with the belt-zone flow pattern were 
then studied incorporating simple heating, to maintain 
the zonal winds, and friction. These objectives were 
evaluated using two values for the Rossby deformation 
length. 
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II. THE MODEL AND FORECAST EQUATIONS 



A. DESCRIPTION 

The atmospheric model for these experiments v;as based 
on the quasi-geostrophic equations. The model was constructed 
on a 3-plane with a central latitude of 45 degrees. Since 
the large-scale motions in the extra-tropical latitudes of 
Jupiter are quasi-geostrophic, the domain of the model was 
restricted to the middle latitudes. Phillips (1963) dis- 
cussed geostrophic motions with approximately uniform 
potential vorticity. These motions have characteristic 
v;avelengths of about 4,000 KM on earth, or the order of 
the Rossby deformation length. He used an expansion in 
the Rossby number, R^ = U/fL, where U is a characteristic 
velocity, f is the Coriolis parameter, and L a character- 
istic length. This expansion showed that the 3-plane 
approximation is justified for geostrophic motions of this 
type. Glerasch, Ingersoll, and Williams (1973) have shown 
that large-scale bands have a characteristic scale of 15,000 
KM in the atmosphere of Jupiter. It was further shown that 

the Rossby number R is much less than one. If the rota- 

o 

tion is sufficiently fast such that R^ is of order one or 
less, Coriolis forces become Important and must be considered 
•in the equations of motion. The additional conditions of 
Rj^/a, where a is the planetary radius, being small, and the 
large-scale waves being greater than or equal to Rj^ are 
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satisfied; therefore, the geostrophlc assumption is valid. 
The geophysical constants for earth and Jupiter are summa- 
rized in the appendix. 

A three-level model was developed so that heating and 
friction could be included and still satisfy the lower 
boundary condition. Since the interior of the planet is 
not knovm, a deep lower layer is used to remove the two 
upper layers from the lower boundary. The solution function 
and its derivatives were forced to vanish at the boundaries. 
The Jovian atmosphere v;as divided into six layers numbered 
1 through 6 as shown in Figure 1. The upper layers, 0-4, 
have thickness AP/2; and the lower layers, 5- and 6, have 
thickness P/2, where P>>AP in this application. The thick- 
ness between levels 3 and 5 is defined to be tt = (AP+P)/2. 

In general these thickness relationships are arbitrary. 

B. METHOD OF SOLUTION 

The quasi-geostrophlc vortlcity equation on a 3-plane is 



If + V-V(c+e^y) - If = 



( 2 . 1 ) 



1 2 

v;here ? = — V (f 



V = ^ kxV(J) 



( 2 . 2 ) 



o 



o 



and the thermodynamic equation is 




(2.3) 
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FIGURE 1 - The Three-Level Model 
Variables are evaluated at the layers shown. 
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dP 

where w = ^ and a is defined as the static stability. The 



quantity t is the radiative time constant and 



iii 

9P 



IS an 



equilibrium temperature which is dependent upon radiation, 
convection, and other small scale parameters. Equation (2.1) 
v;as applied at levels 1, 3, and 5 which yielded 



9? 



3t 



1 ? 

^ . V(?3^+8^y) - , 



o P2 Pq 






( 2 .^) 






^ + v.-v(Co+6.y) - f = 



3t 



’3 o 



o P^-Pp 






3 ’ 



(2.5) 



3^i- Ug W|| p 

- + v^-v(^^+e V) - f = -"2 



3t 



’5 o‘ 



o Pg-Pi, " ^5 > 



(2.6) 



Equation (2.3) applied to levels 2 and 4 yielded 



L- + V .u r ^3~’^i x + „ =1 rri-A _ cM'i ■> 

3t P 3 -P 1 + ^2 ^ ^3"^l ^ 2 2 T ^^3P“^ ^3P^2^ ’ 

(2.7) 




P -P ' 

^5 3 



^ ''V’ ^ “Ii“4 

0 D 



1 [(SLI*!, _ (it), ] . 

T '■^3P ^3P'^V 

( 2 . 8 ) 



Equations (2.7) and (2.8) were solved for Wp and w^j, 
respectively, and substituted with Equation (2.2) into 
Equations (2.^), (2.5) > and (2.6). These operations yielded 
the potential vorticity equations 
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at 



[V^({)^+K2(4>3-4>i>] = -V^-V(v2({)^+fQe^y) - K2V2-V(4>2-<|)3 



■2r. * 



+ Bv (})^ + — [a4>2-(4>3-4>]^) ] = , 



1^ [v^(j)2+Ki,(4.^-«i)2)-K2(<}>3-4'i)] = -V3-v(v^({)2+f^BQy) 






+ BV^ 4>3 '*' ^ 



K. 



- -p[a4>2-(4>3-<{)i)] = C 2 , 



[v2(})5-K5(<j)5-4)3)] = -V5*v(v2<})^+f^6oy) + <^Vi^-v(<p^-<p, 



+ bv^<|.3 - ^[ 3 < j ) 4 -(« J > 3 - 4 ) 3 )] = C3 , 



where 



.2 



^2 *■ .^2 » 



AP a. 



Ki, = 



APua I 



= 6k, 



and 



^6 = 



P.a,^Tr 

D 



= EK, 



lil 



) 

(2.9) 



( 2 . 10 ) 

) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 



A set of three linearly independent eigenvector equa- 
tions were formed to solve the equation set (2.9), (2.10), 
and (2.11) after a method by Charney and Phillips (1953). 
If 



$ 




(2.15) 



then the general form of the potential vorticlty equation 
is 



V 



2 




A K 



2 at 



C 



where 




( 2 . 16 ) 



and A is the appropriate matrix. If 




(2.17) 



is an eigenvector of A and X is its corresponding eigenvalue, 
then 

Ay = X^ , (2.l8) 
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and 



e (2.19) 

Now if 

^ ^ 2^2 ^3)^3 » (2.20) 

and C = + ^ 211.2 ■** '^ 3 li -3 » (2.21) 

are substituted into Equation (2.l6) and Equation (2.l8) is 
used to determine the eigenvalues, then the coefficients of 
each eigenvector can be equated to give 

[(V^+X^K2)H^3 = d^ , (2.22) 

where 1 = 1 , 2 , 3 * 

The following eigenvalues as specified by Equation (2.18) 
were obtained, 

/ 

= 0 

X 2 = - I (2+6+e) + I [(e+6)2+li(i-e)]^ 

X 3 = ’ - ( 2 + 6 +e) - I C(e+6)2+4(l-e)]^ 
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The eigenvector equation set was formed from the poten- 
tial vorticity equations (2.9), (2.10), and (2.11) and 
Equation (2.22). These equations are 



|_ = -V*[e^ykxVH^+V^V%^+V3V^({)3+V^v2|<{.5] 

+ BV^H^ + Q , 



( 2 . 2 ^) 



^ CC7 ^ -V- (g^ykxVH2+V^V^(J)^+e2V3V^4)3+Y2V5V^({)^ 

(2.25) 



C(V^+X3K2)H3] = -V{$^ykxVH3+V^V^(j)^+B3V3V^(j>3+Y3V5V^(J)^ 

- X3K2CV2 ^^^^5 + Q > 

( 2 . 26 ) 

where Q = ^ [ (1-3) (5({)2-8<{)2 ) + (33-ye ) ) 1 • (2.27) 



An assumed solution of the form 

H = H(y,t) + H^cos kx + sin kx , (2.28) 

c s 

where H is the zonal mean solution amplitude, is the 
Fourier cosine eddy amplitude, and H is the Fourier sine 

o 

eddy amplitude, was introduced into Equations (2.2^1), (2.25), 
and ( 2 . 26 ) as suggested by Lorenz (i 960 ). The equations 
vjere quasi-linearized by eliminating all terms of higher 
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wave number than k. This technique permits full interaction 
of the disturbance in x with the mean flow, but neglects 
interaction with the other waves in x. 

Equations (2.2^1), (2.25), and (2.26) were separated 
into components independent of x, coefficients of cos kx, 
and coefficients of sin kx. The resulting equation set 
(2.29) through ( 2 .^ 13 ) constitute the prediction equations 
for this model. The primed quantities represent the partial 
derivative with respect to y. 

The terms Independent of x in (2.2^) give: 



8t 



It 

H 



where 



and 



ft! If! tf f ff f 






r III III X II I II II 

+ - (4>r 4’c: -4>ir 4’cr )~ “ (<Pr <Pc: -4>ir 4>c: )] 
e ^^5c^5s ^5s^5c c ^5s^5c ^5c^5s'^-^ 



+ F + Q 



_ " 

F = BH , 

K 



(2.29) 

(2.30) 

(2.31) 
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The cos kx terms In (2.24) give: 



2 

9 

9t \ 2 

9y 


p Vr p_! 1 n 

_MI p 1 t It 

+ (4>2 +k 4)2<J>2g-4)2'f^3s) 


v/here 


— -ttT p_! 1! 

+ - «!>5 +k <J)5^3^-^3^3^)] + P, + Q, , 

(2.32) 

= B(H”"-2k^H%k^H^) , (2.33) 


and 


K 

-(2.34) 



The sin kx terms in (2.24) give: 



9 .9^ 

9t \ 2 

3y 


- k2)H3 = - 

o 

+ (f3"4'33+k273*3<,-^4'3,) 


v:here 


..11T p_T _l It 

+ F (‘•’5 ®s> 

(2.35) 

11 ft P It il 

Ps = B(Hg -2k^Hg+k^Hg) , (2.36) 


and 


Qg = -f [(l-3)(9<j>2-H2g) + (<53-Ye)0<l>4-94>i|g)] .(2.37) 
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The terms independent of x in (2.25) give: 



It ^jy 2 ^ 2 *^ 2 ^^ 



Ir fit ftl It f II f 



tit til II t II t 

+S2[('t>3c't’33-4'33*3c)-('t3s't’3e-4’3o'<'3s>3 



ttl Itt II t II t 

+Y2t('l’5c*5s-'*'5s'<’5c^''*5s'''5o'%c'*'5s>3 






-’'2<'*’5o'*’3s''*’5s*3o>''^2''’’3o'’'5s"'*’3s*5o^^’ 



+ P + Q 



(2.38) 



The cos kx terms in (2.25) give: 



9 /9 , 2, , XTT _ 

9t 2 
oy 



Ir — Mt p_t _t tt 

^ '*1 h*ls-h*ls> 



— ttt p_t _t tl 

+e2(4>3 <|)33+k 4.34.3g-4)^4.3g) 



_ttt p_t t II 

+Y2(4>5 ^53+k (|)5<f)5^-(^5(|)5^) 



-«t _t 



— t _t 



+ ^2’'2f''*'l'^3s-'*’3'*'ls^''’'2^''’5*3s-'’3'*'5s'^^ ®c’ 

( 2 . 39 ) 
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The sin kx terms in (2.25) give: 



3t ^gy2 ^ ■^^2^2^^s 



(6 f H +((}>; 4), +k'^4>T<}>T -4>-,4>n ) 
f o o c ^Ic ^l^lc 



Mf p — T t »» 

+62(4.3 4.3^+k 4.343^-4.34.3^) 



MT 2 ’ » 

+Y2(4>5 '^5'^5c“'*'5‘^5c^ 



_! t 



t 



+ X2-:2[('t'i'f’3e-'<’3*lo' ^2<'*'5''’3c-*3'‘'5o > ^ • 

(2.40) 



The terms independent of x in (2.26) give 



3 ,3 +7r 

3t ^3y2 



V ttt ftl t! t t! T 

2T” ^^‘^1c'^1s“'^1s‘^1c^“^*^1s‘^1c*'‘^1c'*'1s^ 



ft! ft! ft t It t 

+ 33C(4>3c4>3g-4>3g4>3c)-(4'3s‘^’3c“‘^3c‘*’3s^^ 



m fit tt ! If t 

+Y3C(45^4.53-4.534.5<,)-(4.534.5<,-4.5^4.53)] 



-*‘3''2'^'*3c*ls"*'3s*lc>-^''’ic*3s-*is'*’3o' 



■^3^%o'^3s"*5s*3c>'''^3‘'*’3c*5s"*3s*5c>’) 



+ F + Q , 



(2.41) 
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The cos kx terms in (2.26) give: 



It 'f? ' 

OJ 



k . tn p_f _f tf 

— {3_f* H + ( rfi (f) +k A, d)T -d), (J)t ) 

o o s ^1 ^Is ^l^ls ^l^ls 









.— » — f 



■^^3'^2[(‘J’i^3s'‘^’3‘^’1s^'^Y3(<|)54>3s-4>3(1)5s)]} + ^c > 

(2.42) 



The sin kx terms in (2.26) give; 



|_ (9 -k2+X K2)H^ =. 

9y 






■^®3^*r*3o*’'^'*’3'*'3c‘‘*’3'’’3c' 



+Y3(<f>5 4>5^+k 






f f 



'^^3'^2*- (‘^l‘*’3c ‘^3‘^lc^'^^3^'^5‘^3c"'^3‘^5c^^^ ^s *^s * 

(2.43) 



A method developed by Richtmeyer (1957) was used to 

9H 

solve for the tendency, , at each time step. The equa- 
tions were integrated in time by a marching process. All 
finite differences were centered in time or space. A 
forward time step vias used to start the marching process. 

A finite difference scheme developed by Matsuno (1966) was 
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used every 50 time steps. This is a two step iteration to 
simulate the backward difference method. This scheme was 
applied in this situation to suppress the computational 
mode that would arise from using the centered time differ- 
encing exclusively. A time step of 30 minutes was used. 
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III. RESULTS 



Two experiments were conducted under varied conditions. 
The first experiment was designed to yield the growth rate 
of the most unstable wave in the disturbance fields for 
Rossby deofrmation lengths of 75 KM and 750 KM. This range 
of values for the Rossby deformation length was suggested 
by Gierasch (1973). The upper two levels each had a pressure 
differential of 50 CB, while thickness of the lowest level 
assumed values of 50, 500, and 5000 CB, respectively. The 
y-wavelength of the mean flov/ was 15,000 KM for all cases. 

The x-wavelength varied from 1,000 KM to 15,000 KM for 
Rj^ = 750 KM and from 300 KM to 3000 KM for Rj^ = 75 KM. 

<2 as defined by Equation (2.12) was scaled to be 




The static stability at level 2, 02 > was obtained from 
Equation (3.1) and the static stability at level 4, 
was found by the relation 




where P 2 and are the absolute pressures at levels 2 and 
4, respectively. This ratio was a constant equal to 0.25 
for all experiments. The mean v/lnd was Initialized as 
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a sinusoidal function of y with an amplitude of 20.5 

meters per second in level 1. This value was obtained 
2 

from which is the limiting value for barotrophic 

stability as suggested by Ingersoll and Cuzzi (1969). The 
mean wind was zero in levels 3 and 5. The mean field was 
held constant, which makes the equations linear, throughout 
the first experiment. The disturbance field cos kx contained 
an initial sinusoidal variation with the same y-scale as 
the mean flov; and an amplitude an order of magnitude smaller 
than the mean field. The initial sin kx field was entirely 
described by the cos kx terms. There were no heating or 
friction terms included. 

The forecasts were carried forward until the growth 
rates became exponential. Growth rates of the disturbance 
were calculated by 



1 ^ 
A 9t 



M 




(3.3) 



where M is the exponential growth rate, A is the amplitude 

of the disturbance, and t 2 and t^ are the times at which 

the disturbance amplitudes were A 2 and A^, respectively. 

The maximum growth rate for = 750 KM occurred at an 

x-wavelength of 6,000 KM for all lovrer layer thicknesses 

(Figure 2). For lower layer thicknesses of 50 and 500 CB 

and R, = 75 KM the maximum growth rate was at an x-wavelength 
L» 

of 600 KM. The 5000 CB case showed maximum growth rate at 
an x-wavelength of 500 KM (Figure 3). 



25 . 
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FIGURE 2 - Growth Rate of Y-Dlsturbance Wave for 
Rj^=750 KM 
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FIGURE 3 - Growth Rate of Y-Disturbance Wave for 
Rj^ = 75 KM 
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The phase and amplitude relationships of the disturbance 
field v;ere determined by 

<1) cos kx + ({) sin kx = C cos(kx - 6), 
c s 

where <p and (}> are the geopotentials in the cos kx and 
sin kx fields, respectively. Equation (3.^) was expanded 
to yield the following expressions for the amplitude C 
and the phase angle <5 , 

6 = arctan (~) , (3-5) 

^c 

*^c ‘^s 

C = 2^ = -g.. ^ . (3.6) 

The amplitude and phase relationships for (}>^, and 

4>^ were computed at the maximum growth rate for each case. 
The results are shown in Figures 4 through 15. All ampli- 
tude profiles showed greater relative grov;th in the northern 
wave; however, the effect was more pronounced at = 750 
KM. This was due to the presence of the planetary vorticity 
gradient 3. The relative amplitudes of and 

were the same order of magnitude for the lower level thick- 
ness P = 50 CB (Figures 4 and 7). The amplitude of 
was an order of magnitude smaller than or for 
P 500 CB (Figures 5 and 8). At P = 5000 CB (Figures 6 
and 9) the amplitude of <P^ was negligible compared to 
or 4>2* These characteristics were observed at both values 
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FIGURE 4 - Relative Amplitude of 4>^, ‘f ’5 

as a Function of Y for AP=50 CB, 

Rj^ = 75 KM, and x-wavelength of 6 OO KM. 
Y-Axis has been shifted to the north. 
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FIGURE 5 - Relative Amplitude of (f)^, and 

as a Function of Y for AP = 50 CB, 

P = 500 CB, = 75 KM, and x-wavelength 
of 600 Km. 
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FIGURE 6 - Relative Amplitude of (t)^, and 

as a Function of Y for AP = 50 CB, 
P = 5000 CB, Rj^ = 75 KM, and 

x-wavelength of 600 KM. 
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FIGURE 7 - Relative Amplitude of and 

-L 0 0 

as a Function of Y for AP = 50 CB, 

P = 50 CB, Rj^ = 750 KM, and 
x-wavelength of 6000 KM. 
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-YdO^KM) 

FIGURE 8 - Relative Amplitude of and 

as a Function of Y for AP = 50 CB, 
P = 500 CB, Rj^ = 750 KM, and 
x-wavelength of 6000 KM. 



33 




FIGURE 10 - Phase Angles of c})^, and as a 

Function of Y for AP = 50 CB, 

P = 50 CB, = 750 KM, and 
x-wavelength of 6000 KM. 
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-Y (10^ KM ) 

FIGURE 11 - Phase Angles of (J)^, and as a 
Function of Y for AP = 50 CB, 

P = 500 CB, = 750 KM, and 
x-wavelength of 6000 KM. 
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r^C-URZ 12 - ?h£se Angles of and as a 

Function of Y for A? = pO C3, 

? = 5000 CBj H. = 750 KK, and 
x-wave length of 6OOO KK. 
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FIGURE 13 - Phase Angles of <{>]_» 4 > 2 > a 

Function of Y for AP = 50 CB, 

P = 50 CB, = 75 KM, and 
x-wavelength of 6OO KM. 
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FIGURE - Phase Angles of and <p^ as a 

Function of Y for AP = 50 CB, 

P = 500 CB, = 75 KM, and 
x-wavelength of 600 KM. 
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FIGURE 15 - Phase Angles of (f>^, and <p^ as a 

Function of Y for AP = 50 CB, 

P = 5000 CB, = 75 KM, and 
x-wavelength of 600 KM. 



^0 



of the Rossby deformation length. The phase angles for 
each level showed that all of the waves were barotroplcally 
stable, since the phase lag in y was in the same sense as 
the horizontal wind shear. Also, the vertical tilt was in 
the opposite sense to the vertical shear as is to be ex- 
pected for baroclinlc instability. Note that the y-scale 
of the disturbance is much smaller for = 75 KM compared 
to Rj^ = 750 KM. This result was predicted by Stone (1969) 
who showed that in a slowly varying wind field the y-scale 
of baroclinlc waves V70uld be of the order of the Rossby 
deformation length. It vras observed that the amplitude 
in the disturbance fields for R^^ = 75 KM was contained in 
two narrow regions centered about 2600 KM and 7500 KM. 
Friction and heating were added to the model in the 

second experiment to study the nonlinear effects. The 

^ ^ ^ 

equilibrium temperature profile was obtained from 

the initial value of . This has the effect of main- 

taining the mean wind, since the heating will tend to 
restore the mean temperature field to its initial conditions 
when the disturbance amplitude is sufficient to cause a 
departure. Heating was present in the mean field only. 

The radiative time constant t was 10 Jupiter days. The 

friction term was present in all fields. The friction 

5 5 

coefficient B was equal to 1.0 x 10 and 5.0 x 10 for 
Rossby deformation lengths of 750 KM and 75 KM, respectively. 
Since the grid mesh distance was 150 KM and the x-wavelength 
was 600 KM, the grid may have difficulty resolving the wave; 



therefore, a larger value of the friction coefficient v/as 
needed to maintain stability in the model. The mean field 
was allowed full interaction with the disturbance fields. 

Two cases corresponding to the Rossby deformation lengths 
v;ere performed with a lower layer thickness of 500 CB. 

All other initial conditions were as described in the linear 
experiment . 

The model was Integrated for 104 Jupiter days with a 
Rossby deformation length of 75 KM. At this time the model 
became unstable. The corresponding mean zonal v/ind profile 
and disturbance amplitudes and phase relationships are 
shown in Figures l6 through l8. The maximum wind speeds 
were 25-9 meters per second for the northern easterlies and 
26.3 meters per second for the southern v/esterlles. Zonal 
velocities also developed in layers 3 and 5 - The relative 
minima associated with the maximum wind is the result of 
the amplitude of the disturbance fields. These amplitudes 
are similar to the linear case except the peaks are now 
blmodal and the profile has a broader base, indicating the 
disturbance energy has been distributed over a larger 
region. The phase angles showed barocllnic instability 
in the vicinity of the maximum mean zonal wind in the 
north and the south, since the sense of the vertical wind 
shear was opposite the tilt of the phase. 

Because the 750 KM Rossby deformation length case has 
a much slower disturbance growth rate, the experiment was 
run for 536 Jupiter days. The mean zonal wind profile is 
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FIGURE 16 - Mean Zonal Wind for AP = 50 CB, 
P = 500 CB, = 75 KM, and 
x-wavelength of 600 KM for 
Jupiter Day 104. 
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FIGURE 17 - Relative Amplitude of 0^, (P^, and (p^ 
as a Function of Y for AP = 50 CB, 

P = 500 CB, = 75 KM, and 
x-wavelength of 6 OO KM for 
Jupiter Day 10^1. 




FIGUflE l8 - Phase Angles of and as a 

Function of Y for AP = 50 CB, 

P = 500 CB, Rj^ = 75 KM, and 
x-wavelength of 600 KM for 
Jupiter Day 104. 
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shovm in Figure 22. Relatively small changes have occurred 
compared with the other nonlinear case. The mean wind 
speed is 20.7 meters per second in level 1. Small amplitude 
winds have developed in levels 3 and 5; however, no appre- 
ciable wind is observed in level 5 in the south. Figure 21 
shovjs the amplitude of the most unstable wave in the distur- 
bance field as a function of time. The growth rate is 
exponential with value 

1 9A A T a“6 “1 

^ ^ = 0.21 X 10 sec 

thus the field has not reached quasi-equilibrium. The 
disturbance amplitudes shown in Figure 19 indicate the 
effects of the planetary vorticity gradient 3 have been 
reduced. The northern disturbance has a relatively greater 
amplitude in all three levels. The relative amplitude of 
is an order of magnitude smaller as observed in the 
linear case. 

One inconsistency was present in this experiment which 
has not been resolved. The growth rates observed in the 
linear mode were an order of magnitude larger than the 
nonlinear mode. This was not a physical result, but an 
error in the program. 
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-Y (10^ KM) 

FIGURE 19 - Relative Amplitude of (j)^, (J)^, and 
as a Function of Y for AP = 50 CB, 
P = 500 CB, Rj^ = 750 KM, and 
x-wavelength of 6000 KM for 
Jupiter Day 536. 
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FIGURE 22 - Mean Zonal Wind for AP = 50 CB, 
P = 500 CB, = 75 KM, and 
x-wavelength of 600 KM for 
Jupiter Day 536. 
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IV. CONCLUSIONS 



A three-level quasl-geostrophic atmospheric model has 
been examined on a 3-plane with a single wave In the x- 
dlrectlon and a deep lower layer of varying thickness. The 
disturbance structure In y was Investigated to determine 
the linear Interactions vjlth a constant mean wind and for 
nonlinear Interactions with a varying mean wind. Simple 
heating to maintain the mean wind and friction were Included 
In the nonlinear case. Two values of the Rossby deformation 
length were Included In the experiments. 

The linear experiment showed that the growth rates In 
the disturbance field are largely Independent of the lower 
layer depth (Figure 2). These growth rates are a function 
of the Rossby deformation length and the x-wavelength . The 
waves v;ere barotroplcally damped In this case, since the 
relative vortlclty gradient was less than or equal to the 
planetary vortlclty gradient. This was obvious since the 
phase tilt v;as In the same sense as the horizontal wind 
shear. 

The heating function v/as not Important In the nonlinear 
experiment, since the model never reached a long term steady 
state. The disturbance fields were barotroplcally damped. 

In the cases studied the zonal disturbances are barocllnlcally 
unstable. If the actual wind structure on Jupiter Is baro- 
cllnlcally stable as suggested by Ingersoll and Cuzzl (1969), 
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then this could be accomplished by distributing the mean 
change through great depth. This would give a very small 
shear which could be stabilized by B. The growth rates were 
expected to decrease with increasing lower layer depths, 
but in fact they did not. Thus it may be that the lower 
layer was not properly formulated by this model. 

Further investigation of the disturbance growth rates 
is required to resolve the discrepancy between the linear 
and nonlinear cases. It is also recommended that the grid 
mesh be reduced to more realistically model the 600 KM 
x-v;avelength. This should allow use of a smaller friction 
coefficient . 
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APPENDIX 



GEOPHYSICAL AND ASTRONOMICAL CONSTANTS OF JUPITER AND EARTH 



PARAMETER 


JUPITER 


EARTH 


Perihelion 
Distance, A.U. g 

(1 A.U. = 1-.496 X 10 KM) 


^.951 


0.9833 


Amphelion 
Distance, A.U. 


5.^55 


1.0167 


Siderial 
Period, (day) 


10 , 759.20 


365.256 


Inclination of 
orbit to ecliptic 


3°04 ' 


230271 


Mass, (GM) 


1.899 X 10^° 


5.976 X 10 


Density, (GM/CM^) 


1.33 


5.517 


Mean Radius, (KM) 


69,828 


6,371 


Equatorial 

Gravitational P 

Acceleration, (CM/SEC ) 


2587 


979 


Albedo 


0.73 


0,36 


Orbital 

Eccentricity 


0.0^183 


0,01627 


Mean Temperature 


130° K 


330° K 


Solar 

Constant, (Ly,MIN) 


-0.07 


-2 


fi, (SEC"^) 


17.6 X 10 "^ 


7.3 X 10 "^ 


Observed Mean 
Vfind, (M/SEC) 


50 


15 


Rossby Number, R^ 


O.OOill 


0.032 


Rossby 

deformation Length, (KM) 


" 765, 78.5 


1,000 
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PARAMETER 



JUPITER 



EARTH 



Coriolis Parameter, (^5°) 

Planetary Vorticity_-, 
Gradient , ), (SEC 

R^/a 



2.5 X 10 ^ SEC 

3.6 X 10"^^ 



1.0 X 10 ^ SEC 

1.6 X 10 



0.05 



0.63 
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